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Abstract 

We construct models of leptons based on 5*4 family symmetry combined with a gen- 
eralised CP symmetry Hcp- We show how the flavon potential can spontaneously 
break the symmetry Si xi Hcp down to Z2 x H^p in the neutrino sector, where 
the choice of preserved CP symmetry H^p is controlled by free (real) parameters 
in the flavon potential. We propose two realistic models of this kind, one at the 
effective level and one at the renormalisable level. Both models predict trimaximal 
lepton mixing with CP being either fully preserved or maximally broken, with the 
intermediate possibility forbidden by the structure of the models. 
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1 Introduction 



Recent hints [1-3] for a nonzero reactor mixing angle ^13, and its subsequent measurement 
by the Daya Bay [4], RENO [5] and Double Chooz [6] Collaborations have revealed an 
unexpectedly large mixing angle of about 9° to an accuracy comparable to the other 
two lepton mixing angles, i.e. the solar and atmospheric angles [7]. This observation 
has had a major impact on models of leptons based on discrete family symmetry (for a 
recent review see e.g. [8]). On the experimental side, this remarkable discovery opens up 
the possibility of unravelling (some of) the remaining unknown parameters of the lepton 
sector, the neutrino mass ordering as well as the leptonic CP violating phases,^ which is 
the primary goal of next generation neutrino oscillation experiments. 

Concerning the Dirac CP violating phase 5c p-, the latest global analyses of the available 
neutrino oscillation data do not allow to pin down a preferred value at the So" level [9-11]. 
Yet, as we move into the age of precision measurement of the lepton mixing parameters, 
this state of affairs is bound to change. Inspired by the successes of discrete family sym- 
metries in predicting the lepton mixing angles, it is not unreasonable to believe that the 
symmetry approach can also lead to predictions for 5c p- To find a realisation of this idea, 
it appears natural to combine a family symmetry and generalised CP symmetries such 
that both the mixing angles and the CP violating phases are predicted simultaneously. 
In this setup, the breaking of the family symmetry would give rise to particular mixing 
angles, while the breaking of the generalised CP symmetries leads to particular values of 
the CP phases. 

The concept of generalised CP transformations has been around for decades, and it was 
realised early on that there exists a one-to-one correspondence between the generalised 
CP transformations and the automorphism group of an imposed family symmetry [12-16]. 
Recently, the consistent generalised CP transformations (i.e. the automorphism group) 
of the discrete family symmetry groups of order smaller than 31, with irreducible triplet 
representations, have been discussed in [17]. 

Despite these fundamental studies, there has been only little work on constructing 
explicit models which realise the idea of combining family and generalised CP symmetries 
[18-20]. In Ref. [19], an 5*4 model with imposed generalised CP symmetry is presented 
which, adopting the type II seesaw mechanism, predicts a normal neutrino mass ordering, 
a maximal atmospheric angle and a maximal Dirac CP violating phase, 5c p = ±7r/2. 
Based on pure symmetry arguments, a systematic classification of an S4, family symmetry 
combined with generalised CP symmetries is performed in [20]. Although the resulting 
phenomenology was analysed in great detail, a dynamical model which implemented these 
ideas was not given in [20]. 

In this paper we shall address the question of breaking a generalised CP symme- 
try spontaneously in the framework of models based on 6*4 family symmetry. We shall 
construct models of leptons imposing both 6*4 family symmetry and the corresponding 
generalised CP symmetry Hcp- We shall then show how the flavon potential can sponta- 

^Thc leptonic CP phases include one Dirac phase 5cp and two Majorana phases 0121, a^i- The former 
can be measured in neutrino oscillation experiments, the latter are relevant for ncutrinoless double beta 
decay and more difhcult to be measured. 
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neously break the symmetry 5*4 xi Hcp down to Z2 x H^p in the neutrino sector, which was 
simply assumed to happen in [20]. In our models, the choice of preserved CP symmetry 
H^p in the neutrino sector is controlled by free (real) parameters in the flavon potential. 
We propose two realistic models of this kind, one at the effective level and one at the 
renormalisable level. Both models predict trimaximal lepton mixing with CP being either 
fully preserved or maximally broken, with the intermediate possibility forbidden by the 
structure of the models. 

The layout of the rest of this paper is as follows. In section 2 we draw together some 
of the basic observations about generalised CP symmetry when combined with a discrete 
family symmetry that are rapidly becoming to be appreciated in the literature. In section 3 
we focus on the case of 5*4 family symmetry, which was also analysed extensively in [20]. 
However, in our case, we will work in a different basis, the diagonal charged lepton basis, 
where the representation matrix for the order three generators T is diagonal, and we verify 
explicitly that the same physical results emerge, as expected. This discussion also sets 
the notation and allows us to be rather brief when discussing the physical implications 
of the models which follow. In section 4 we present our first example of an 5*4 x Hcp 
model at the effective level (i.e. involving non-renormalisable operators) where the fiavon 
potential leads to exactly the kinds of possibilities for CP violation discussed in section 3. 
In section 5 we present a superior renormalisable 5*4 x Hcp model which has the added 
advantage of explaining why the reactor angle is smaller than the other lepton mixing 
angles, and also leads to similar options for CP violation. Section 6 concludes the main 
body of the paper, followed by two appendices on more technical aspects. 



2 A consistent definition of generalised CP transfor- 
mations 

In this section, we start by briefly reviewing for completeness how CP transformations are 
generally defined in a consistent way without an existing family symmetry. This is begun 
by noting that in the extension of the Standard Model (SM) with Majorana masses for 
the left-handed neutrinos, the mass terms and the charged current interactions for the 
lepton fields read (below the electroweak symmetry breaking scale), 

£ = -hmiln - ^ulCm.UL + Li^^lW- + h.c. (2.1) 

= (e, /i, r)^^ stands for the SM charged lepton fields, = {Uf., z/^, i^r)I are the left- 
handed neutrino fields, C is the charge conjugation matrix, mi^u are complex matrices, 
and is symmetric. The charged current interactions are invariant under the so-called 
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generalised CP transformation defined by^ 



Mx) ''4iXli'CuI{x'), (2.2) 

ln{x) ^2X^7°C5(x'), 

wfiere and are unitary matrices acting on generation space and x' = (t, — x). 
Notice that the "canonical" /traditional CP transformation can be recovered from the 
above equations by setting both X^ and Xji to the identity matrix. Furthermore, note 
that the lepton fields II and have to transform in the same way due to the presence of 
the left-handed charged-current interactions. Then, the Lagrangian of Eq. (2.1) conserves 
CP if and only if the lepton mass matrices mi and rrij, satisfy the following relations, 

XlmiXR = ml Xlm,XL=ml. (2.3) 

Let us now consider a theory that is invariant under both a generalised CP symmetry 
and a family symmetry Gf, which is the main focus of this paper. We include this 
discussion in order to set the notation and to make the rest of the paper accessible. 
However, since all the information in this section is already in the literature, the discussion 
below is necessarily brief and we refer interested readers to e.g. [17,20] for more details. 
Let us, then, consider a field in a generic irreducible representation r of G/ which 
transforms under the action of as 

G 

^{x) PrigMx), geGf , (2.4) 

where pr{g) denotes the representation matrix for the element g in the irreducible repre- 
sentation r. The mapping of (p under a generalised CP transformation is given by [12-14, 
17,20,21]: 

ip{x)^X,ip*{x') , (2.5) 

where Xj. is a unitary matrix in order to keep the kinetic term invariant, and the obvious 
action of CP on the spinor indices has been omitted for the case of ip being a spinor. 

Requiring a generalised CP symmetry Hcp in the context of a flavour model based 
on some family symmetry Gf, restricts the allowed choices for X^ considerably [17,20]. 
If we first perform a CP transformation, followed by a family symmetry transformation, 
and subsequently an inverse CP transformation we obtain 

^ X,v*{x') % X,p:{g)^\x') ""-^ X,p:{g)X;'^{x) . (2.6) 

As the theory should be invariant under this sequence of transformations, the resulting net 
transformation must be equivalent to a family symmetry transformation Pr{g') of some 
group element g' 

X,p:{g)X;' = p,{g'), g,g'eGf. (2.7) 

factor of i is included, differing from the definition of Ref. [21], in order to eliminate an extraneous 
"— " sign from the second equation of Eq. (2.3). 
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This equation defines the so-called consistency equation which must be satisfied in 
order to implement both generalised CP and family symmetries simultaneously. It is 
important to note that Eq. (2.7) must hold for all representations r simultaneously, i.e. 
the elements g and g' must be the same for all irreducible representations. Furthermore, 
Eq. (2.7) implies that the generalised CP transformation Xr maps the group element g 
onto g', and this mapping preserves the family symmetry's group structure.^ For faithful 
representations r where the function pj. maps each element oi Gf into a distinct matrix, 
the consistency equation will define a unique mapping of the abstract group Gf to itself.^ 
It is also clear from Eq. (2.7) that g and g' must be of the same order. 

It is important to note that, if Xr is a solution to Eq. (2.7), then not only is e^^Xr a 
solution (with arbitrary phase 6), but also 

pr{h)Xr, with heGf, (2.8) 

are solutions as well.^ Therefore the consistency equation can only determine the possible 
form of the CP transformation Xj. up to an overall arbitrary phase and a. Gf transforma- 
tion Pr{h) for a given irreducible representation r. This statement applies to any family 
symmetry group Gf regardless of it being discrete or continuous. For discrete family 
symmetries, it is sufficient to impose the consistency equation Eq. (2.7) on the group's 
generators. 



3 84^1 Hcp scenarios with Z2 x H^p preserved in 
neutrino sector and Zs in the charged lepton sector 

In the remainder of the paper, we will investigate the generalised CP transformations 
Hcp consistent with 5*4 as a family symmetry group, 

^4 X Hcp. (3.1) 

In other words we shall be interested in theories which (at high energies above the sym- 
metry breaking scale) respects 5*4 x Hcp where the fields transform as 

V9(x) ^ pr{g)ip{x), g e S4 , (3.2) 

¥^{x) ^ Xr^*(x'), (3.3) 

^X,p*{g,g2)X-^ = X,p*jg,)p*{g2)X-^ = XrP*{gi)X-'X,p*{g2)X-^ = p,{g[)p,{g'^) ^ p,{g[g'^), 
where we denote XryO*((?i)Xj. ^ = Pr{g'i)- Therefore the CP transformation X^- is a homomorphism of the 
family symmetry group Gf. 

"^Thc consistency equation defines an automorphism for tlie group Gf. See Ref. [17] for a more formal 
treatment. 

^To see this note that pr{h)X„pl{g){py{h)Xy)-'^ = py{h)Xypl{g)X-^ p-^{h) = p„{h)py{g')py{h-^) = 
Pr{hg'h~^). Therefore the generalised CP transformation pr{h)Xr maps the group element g into hg'h^^ 
which belongs to the conjugacy class of g'. It is equivalent to a CP transformation X^ plus an inner 
automorphism of the family symmetry group Gf [17]. 
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where the elements of Hqp are denoted by X^. Since 6*4 is our primary focus, in Ap- 
pendix B, we prove exphcitly the statement in [17] that the most general CP transfor- 
mation Hcp consistent with the 5*4 flavour group is simply given by being equal to 
the identity (up to an inner automorphism). Therefore the most general CP transforma- 
tion Hcp consistent with 5*4 family symmetry is of the same form as the family group 
transformation itself, 

Xr = Pr(/i) , he Si, (3.4) 

where h can be any of the 24 group elements of 5*4. 

Following [20], we further assume in this section that the underlying combined sym- 
metry group 5*4 XI Hcp is broken to G^p in the neutrino sector and G^ in the charged 
lepton sector^ 

G^cp = ^2 X if^p and = Z3 . (3.5) 

The main purpose of this paper (discussed in the following sections) is to show how the 
original symmetry 6*4 xi Hcp can be spontaneously broken to Z2 x H'^p in the neutrino 
sector, which goes beyond the analysis in [20]. However, in this section, we begin by 
following in the footsteps of [20] but using a different 6*4 basis, namely that in which the 
charged leptons are diagonal and all the mixing arises from the neutrino sector. In fact we 
regard it as necessary to repeat the symmetry analysis in the new basis in order to verify 
explicitly that the two bases give equivalent results, thereby placing our later results on a 
sound footing. It is also useful to set the notation and to allow us to arrive at the physical 
results very efficiently in the next sections. Therefore, although the results in this section 
may seem like repetition, we regard it as both necessary and useful to proceed by first 
following the analysis in [20] in a different basis. 

The basis in which we work in this paper is motivated by the trimaximal 6*4 model 
of [22] (originally proposed without Hcp), where the resulting symmetry breaking led to 
a preserved Zf in the neutrino sector generated by S, as well as a diagonal charged lepton 
mass matrix which arises thanks to the chosen 6*4 basis with diagonal complex T generator, 
as discussed in Appendix B. The Z2 symmetry of the neutrino sector will automatically 
produce a neutrino mass matrix in which the second column of the corresponding mixing 
matrix is proportional to (1, 1, 1)"^. In the next section we shall discuss models which 
are inspired by this model but include Hcp and are hence based on 5*4 x Hcp which is 
spontaneously broken to 

Gcp = Z2 X H^p , (3.6) 

in the neutrino sector. For G^p to be a well-defined symmetry, the consistency relation 
of Eq. (2.7) needs to be satisfied for the residual symmetry group . In other words the 
elements of H^p must satisfy, 

Xrp:{S)X-' = pr{S). (3.7) 

One can arrive at the above restricted consistency equation by recalling the general con- 
sistency condition of Eq. (2.7), with g,g' G , and realising that g' must be S, as it is 

^In the models constructed in sections 4 and 5, the Z3 subgroup of 5*4 will be broken. However, this 
happens in such a way that the resulting charged lepton mass matrix remains diagonal. 
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the only element which (trivially) has the same order a.s g = S. Armed with this in- 
formation, recall from section 2 that a faithful representation of the group Gf must be 
used to uniquely determine the mapping of the consistency equation (then the unfaithful 
representations will follow). Hence, we proceed by considering faithful representations 
of 5*4, i.e. the three dimensional irreducible representations. 

For the triplets, one can easily show that for H^p there are only eight choices for the 
Xj. in Eq. (3.7) which are acceptable: 

Xr = Prih) , (3.8) 

with 

he{l , S , TST^ , T^ST , U , SU , TST'^U , T^STU} . (3.9) 

This may be compared to the 24 choices for Xr in Eq. (3.4) corresponding to the original 
Hcp before it is broken (in the charged lepton sector). Note that the condition in Eq. (3.7) 
is automatically satisfied for the doublet and the singlet representations (the unfaithful 
representations), as is implied by the consistency equation. Comparing the results of 
Eq. (3.9) to the six solutions of the analysis given in [20], we observe that our third and 
fourth solutions of Eq. (3.9) are missing in [20]. This is because the authors of [20] have 
made the assumption that a generalised CP transformation must be symmetric. However, 
it is not necessary to demand symmetric Xj. if there is a family symmetry, so we keep all 
eight solutions in the following. 

Let us assume that the three generations of lepton doublets reside in the 3-dimensional 
irreducible representation.'' Then, it is possible to construct the light neutrino mass 
matrix by demanding that it respects both the residual flavour symmetry and the 
generalised CP symmetry H^p-. 

pl{S)m^Ps{S) = my, (3.10) 
Xlm,Xs = ml. (3.11) 

The most general matrix which satisfies Eq. (3.10) has the form 

/ 2 -1 -1\ /I 0\ /O 1 1\ / 1 -1\ 
my = a\-l 2 -1 +/3 1 +7 1 1 +e 1 -1 , (3.12) 
\-l -12/ \0 1 0/ \1 1/ \-l 1 / 

where a, /3, 7, and e are complex parameters. The CP symmetry of Eq. (3.11) will further 
constrain them to be either real or imaginary. The eight possibilities listed in Eq. (3.9) 
are classified in Table 1 according to their four distinct predictions for a, /3, 7, and e being 
real or imaginary. 

In order to diagonalise m^ of Eq. (3.12), we perform a tri-bimaximal transformation 
Utb, 

/ 3a + /3 - 7 -VSe \ 

K = UlB^yUTB=[ /3 + 27 , (3.13) 

V -\/3e 3a-/3 + 7 / 

^The results would be the same if the lepton doublets were assigned to 3', since the representation 3' 
differs from 3 only in the overall sign of the generator U . 
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^CP 


a 


Q 
P 


7 


e 


(i) 


Pr(l), PriS) 


real 


real 


real 


real 


(ii) 


PriTST''), priT^ST) 


imaginary 


real 


real 


real 


(iii) 




real 


real 


real 


imaginary 


(iv) 


p^iTST^U), p^iT^STU) 


imaginary 


real 


real 


imaginary 



Table 1: The generalised CP transformations consistent with a residual Z2 symmetry in the neutrino 
sector and their corresponding constraints on a, /3, 7 and e parameters in Eq. (3.12). 



with 



U' 



TB 



]_ j_ ]_ 

\ L J- J- I 

\ V6 V3 V2 / 



(3.14) 



followed by a (1,3) rotation Ul, 

U'Jm'^Ul = diag(mi,m2,m3). 

Then the PMNS matrix is 

UpMNs = UtbUI , 



(3.15) 



(3.16) 



where Upmns is only determined up to the permutations of row and column, since the 
order of both the charged lepton and neutrino masses cannot be predicted in the present 
symmetry guided approach, see also [20]. Finally, we shall work in the PDG convention [7], 



UpMNS = V^diag(l,e* 2 ,e 



"31 . 



(3.17) 



with 



V 



C12C13 S12C13 

-S12C23 - Ci2S23Sl3e*''^'P C12C23 - Si2S23Sl3e*'^^-P 



Sue 
S23C13 
C23C13 



iScp 



(3.18) 



■S12S23 ~ Ci2C23Si3e — C12S23 — Si2C23Si3e 

We are now prepared to consider each of the four cases of Table 1 in turn. 

(i) H'^p = pr(l), Pr{S) 

In this case, the unitary matrix is 

cos 9 sin 9 
U',.= \ 10 
— sin 9 cos 9 



with 



tan 29 



\/3 e 



(3.19) 



(3.20) 
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The transformation P in Eq. (3.19) is a unitary diagonal matrix with entries ±1 or 
±2, which renders the hght neutrino masses mi 2,3 positive. In the following cases, 
we shall omit this trivial factor. Given the diagonal charged lepton mass matrix, 
the leptonic mixing matrix is 



UpMNs — UtbUI 



I COS 61 



I sin 6' 



V-76^°s^ + 75^i^^ 73 -TS^i^^ ' 



V2 



COS 



As a result, the lepton mixing angles and CP phases take the values 



P . 
(3.21) 



sin^^ 



13 



sina2i = sina3i = sm.5cp = , 
I sin^ e, sin^ = iT^' ^^^^ ^23 = 



2+cos 26* 



(3.22) 

(3.23) 



Notice that there is no CP violation in this case, but it is otherwise a viable scenario. 
Finally, the light neutrino masses are 



mi 



3a - 



7712 = I/3 + 27I 



m3 



3a + V(7 - + 3e^ 



(3.24) 



(ii) H^cp = Pr{TST'),pr{T'ST) 

This case is not discussed in [20] because the authors require that the CP transfor- 
mation be both unitary and symmetric. However, it turns out that for this case 

-9a2 + (/3 - 7)2 + g 

77i',7n'J = \ (/3 + 27)' 

-9a^ + {(3- 7)2 + 3e2 

which implies that |mip = |m3p. Therefore, the resulting light neutrino mass 
spectrum is (partially) degenerate, and the PMNS matrix cannot be determined 
uniquely. Phenomenologically, this case is clearly not viable. 

(iii) H-p = Pr{U), p,{SU) 

The unitary transformation Ul is of the form 

cos 9 sin 9 
K=\ 1 I , (3.25) 



with the angle 9 



-isin6' icos9 



7€ 

tan29 = ^^. (3.26) 
v3 a 
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The PMNS matrix takes the form 

/ JlcosO 



UpMNs — UtbUI 



\ 



..igcos^-^sin^ 



1 








1 


1 












V6 



I sin 9 



-4^sm6 - ^ cos 6 



(3.27) 



V2 



COS 9 j 



Therefore the lepton mixing parameters are 

0, sinasi = 0, 



sm a2i 
1 



I sin b^ 



CP 



sin^ 6*13 = I sin^ 6', sin^ 6*12 



1 



2+cos IB ' 



sin^ 6^23 



(3.28) 
(3.29) 



Notice that we have maximal Dirac CP violation ^cp = if- However, the sign 
cannot be fixed uniquely because it depends on the value of the angle 6 and the 
order of the light neutrino masses. This case provides a maximal atmospheric mixing 
angle, with the light neutrino masses given by 



mi 



/3 - 7 - V9a^ - 3e2 



1712 = I/3 + 27I 



3e2 



(3.30) 



Therefore this case is viable as well. 

(iv) H^p = priTST^U), priT^STU) 

In this case, the unitary matrix Ul is given by 



U' 



1 

7! 



--) 



-e^(^+f) -( 
1 



with 



tan26' 

The resulting PMNS matrix reads 



-/3 + 7 
3ia 



(3.31) 

(3.32) 

\ 

' I 
(3.33) 

Agreement with the present experimental data can be achieved if we permute the 
rows of the above U pmns by [20] 



U PMNS — UtbUI 



(e+f) 



<e-f) 



1 

1 



1 



73^ 



1 1 e<^-f ) + ie^<^+f ) 



7-rlsi 
PMNS 




PMNS, or f/ 



2nd 

PMNS 




PMNS 



(3.34) 



to yield two phenomenologically viable PMNS matrices, U]^lj^g and U]l^^j^g, so 
named because they differ only in their prediction of the atmospheric mixing angle. 
^¥mns predicts the atmospheric mixing angle to be in the first octant and W^^j^g 
in the second. The three lepton mixing angles are then determined to be ^ 



sin^^ 



13 



sm ^23 



1(2 

6 



V3cos2^ 



sin^^ 



or 



4 + 73 cos 26 
Both solutions share the phase predictions 

V3 + 2cos2^ 



12 



^•2 Q2nd 
sm 6^23 



v/3 cos2e 
2 



sina2i 



2 + 73 cos 2^ 



|sin a^^\ 



4 + 73 cos2^ 



4^3 sin 2^ 



(3.35) 



5 - 3 cos 4^ 



I sin (5, 



CP 



^4-2^3 cos 2^ (4 + 73 cos 2^) sin 26 



5 - 3 cos 4^ 



(3.36) 



where a'-^i = a^i — 25c p-, the parameter a'2,1 has been redefined to include 5c p which 
is useful in the context of neutrinoless double beta decay [21]. Lastly, we find the 
light neutrino masses for this case to be 



nil 



a' 



ni2 = I/3 + 27I 



73ie+ 7(/3-7)2 -9a2 



(3.37) 



Therefore this case is also viable and involves non-trivial CP violating phases. 



4 An effective 5*4 >^ Hcp model 



Field 


L 








r= 


Hu,d 




V 















e 




3 


3 


1 


V 


1 


1 


3 


2 


3' 


1 


2 


3' 


1 


3' 


1 


z. 


1 


1 


i 


-1 


—i 


1 


i 


i 


1 


1 


1 


-1 


-1 


1 


1 






u 






u 


1 


1 


1 


u 


CO 


cu 


1 


1 


UJ 


UJ 


U{1)r 


1 


1 


1 


1 


1 




















2 


2 


2 


2 



Table 2: The particle content and the transformation properties under the family symmetry S^x Zj^x 
and U{1)r. 



In this section, we present a realistic effective model of leptons based on 

^4 X Hcp , (4.1) 

^The following results for the mixing angles and CP phases match those of [20] after making the 
replacement 6 7r/4 — 6. 
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supplemented by the extra symmetries 

Z4 X Z3 X Uil)R (4.2) 

in order to control the allowed operators in the model. 

The three generations of left-handed (LH) lepton doublets L and the three generations 
of right-handed (RH) neutrinos A^^ are both unified into the 6*4 triplet representation 3, 
while the RH charged leptons are assigned to be the singlet representations 1 or 1'. The 
complete list of lepton, Higgs, flavon and driving fields as well as their transformation 
properties under the family symmetry are listed in Table 2. Notice that the Z3 symmetry 
is used to separate the flavons entering the charged lepton sector at leading order (LO) 
from those of the LO neutrino sector. They are further distinguished by an auxiliary Z4 
symmetry. This Z4 symmetry is also helpful in achieving the charged lepton mass hierar- 
chies. As discussed in the previous section, the 5*4 family symmetry will be spontaneously 
broken to generate trimaximal mixing. 

Furthermore, we will impose a generalised CP symmetry Hcp consistent with 5*4 as 
discussed in the previous section. In addition to forcing the coupling constants in the 
superpotential to be real, the generalised CP symmetry will also be spontaneously broken 
in the charged lepton sector. In the neutrino sector, a restricted generalised CP symmetry 
remains, 

G^P - Z| X H^cp , (4.3) 

where this residual symmetry results in predictions for Sep and in other relations as 
discussed in section 3. However, unlike in the previous section, the dynamics of the 
spontaneous breaking of S4 x Hcp will be discussed by studying the vacuum alignment 
of a suitable flavon potential to which we now turn. 

4.1 Vacuum alignment 

The general driving superpotential invariant under the symmetry of the model is 

+ /l {fs if Sf 3)3') I + /2 (^S (te)3')l + /s (^SV'S)!^ 

+he {fsfs)^ + h^e (00)i + h^ee , (4.4) 

where (. . .)r denotes the contraction of the indices to the representation r. If we 
require the theory to be invariant under the generalised CP transformation defined above, 
then all the couplings gi, and hi are real parameters. In the SUSY limit, the vacuum 
configuration is determined by the vanishing of the derivative of the driving superpotential 
Wd with respect to each component of the driving fields. The vacuum in the charged lepton 
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sector is determined by 

dWd _ r, , 2 



dwd _ ^ I 2 
dwd 



93(^1^ + Sv^TaV^Ta) + '294VlV2 = . (4.5) 



This set of equations admits two solutions, the first one is given by 

(<^t) = j 1 Ut, = ( 1 ) ^^^^ '"t = -'^^1 , (4-6) 



and the second solution reads 





(ipr) = yi jvT, = ( 1 ) = ^Vr, . (4.7) 

Here, we choose the second solution.^ The phase of can be absorbed into the lepton 
fields, therefore we can take to be real, then the VEV vt is real as well. The equations 
determining the vacuum alignment in the neutrino sector are 

2/l((/5|j - ^32^3^) + f2{(pl^S2 + (l>2^Sz) + MVS^ = , 



dw, 



dwd 



-- 2/i((^|2 ~ VSiVSs) + /2(01<^5i + <P2VS2) + fs^VSs = 
= 2/i(</?|3 - ips,VS2) + f2{(pl^S3 + (t>2^Sr) + M^S2 = 



= h{^l^+2^s2Vss) + 2h2<Pi<P2 + he = 0. (4.8) 

Disregarding the ambiguity caused by the S4 family symmetry transformations, we find 
the solution 

iVs) =\l\vs, {^)=^ll], {0 = n, (4.9) 

where the VEVs obey the relation 

v,+V2 = -f3u/f2 = -Bu , (4.10) 

3hivl + 2/i2t'if2 = -hsu"^ . (4.11) 



^Notice that the first solution can be eliminated by adding another driving field with the same flavour 
numbers as [23]. 
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The VEVs ^5,^1 and V2 are undetermined. In order to solve this problem, we introduce a 
second driving fields with identical quantum numbers as [23] . Consequently we obtain 
two F-term conditions which are identical in their structures but involve independent 
coupling constants, hi for one driving field and h'^ for the other, i.e. 

3h[vl + 2h2ViV2 = -h'^u^ . (4.12) 

As both conditions, Eq. (4.11) and Eq. (4.12), must be satisfied, one can find a unique 
solution for the VEV of the ^ps fiavon 

and 

Since it is always possible to absorb the phase of u by a redefinition of the matter fields, 
we can take u to be real without loss of generality, then vs is either real or purely imagi- 
nary. The vacuum of the doublet </> can be obtained by solving Eq. (4.10) together with 
Eq. (4.14). 

• For fi2 _ 4C> we find 



vi = \[-B± ^B^ - AC ] u 
V2 = \ [-Bt VB^ - 4C ] u 

from which we see that f 1 and f 2 are both real. 
For 52 - 4C < we obtain 



(4.15) 



(4.16) 



vi = \ [-B ± jy /AC - B^ ] u , 
V2 = \ [-B T W^C - 52 ] u , 

which entails Vi = V2 in this case. 
4.2 The lepton masses and mixings 

The effective superpotential for the charged lepton masses which is allowed by the sym- 
metries is given as 

W, = ^(L^^)^if,r^+||(L(^T^T)30v^<i/"' + 7^(^(^^T)30v^<i/i'^ 

i 

where O stands for 

O = {cp^T, wl, rfipr] , (4.18) 
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and all possible 5*4 contractions are to be considered. It is interesting to note that the tau 
mass is suppressed by 1/A, while the muon and electron masses appear at order \/ A? and 
respectively. The mass hierarchies among the charged lepton are thus produced in 
a natural way. Inserting the vacuum alignment of Eq. (4.7), we find a diagonal charged 
lepton mass matrix with 



Vt 



2i/mi (y) -yi^2^^ ^d, me = ye(^) Vd, (4.19) 



where ye is the result five different contributions corresponding to the in Eq. (4.17). 
Obviously the VEVs of the flavons ipx and t] are responsible for the spontaneous breaking 
of both the family symmetry and the generalised CP symmetry. It is straightforward to 
check that 5*4 is broken completely in the charged lepton sector. In addition, only one 
CP symmetry = Pr(l) out of the 24 consistent CP transformations is preserved by the 
flavons ifiT and r], given the previously mentioned fact that both vt and can be chosen 
to be real. As a consequence, the residual CP symmetry in the charged lepton sector is 
H^p = Pr(l)- Now we turn to the neutrino sector, the effective superpotential is of the 
form 

= y {LN\ + y, {{N^N^)^, ^s)^ + 2/2 (N^N^), ^ + ys {{N'N^), , (4.20) 

where all the couplings yi are real due to the general CP invariance. The Dirac neutrino 
mass matrix takes a trivial form 



mo = 1 \yv^. (4.21) 




Given the vacuum conflguration of Eq. (4.9), it is straightforward to derive the RH neu- 
trino mass matrix 




mM = a\ -1 2 -l +6 00 1|+c|0 10|+(i|100|, (4.22) 
\ -1 -1 2 / \ 1 C 

where we have introduced the parameters 

a = yivs, b = y2U, c = ygWa, d = y^vi . (4.23) 

The light neutrino mass matrix is then obtained from the seesaw formula, yielding 

= — m/jm^/m^ 

/I 0\ /O 1 1\ / 1 -1\ 
+ /3 1+7 ll0 +e l -1 . (4.24) 

\0 1 0/ \1 1/ \-l 1 / 

This matrix is of the same form as the neutrino mass matrix in Eq. (3.12) which is the 
most general matrix invariant under the action of the 5*4 element S. The parameters a, 
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a 




7 


e 


(i) 


aeR, beR, ceR, deR 


real 


real 


real 


real 


(ii) 


aeC, 6eM, ceM, deM with a = -a* 


imaginary 


real 


real 


real 


(iii) 


aeR,beR,ceC,deC with c^d* 


real 


real 


real 


imaginary 


(iv) 


aeC,beR, ceC,deC with a = -a* and c = d* 


imaginary 


real 


real 


imaginary 



Table 3: The four possible cases of the effective model which are controlled by the input parameters of 
the flavon potential. 



/3, 7 and e are related to a, b, c and d by 

a 

a = 



f3 
7 



-9a2 + 62 + (c + rf)2 - 3cd - b{c + d) ' 
Sa^ - 6^ + 

ih + c + d) [-9a2 + 62 + (c + df - 3cd - b{c + d)] ' 

6a2 - (c + (i)2 + 2cd + 6(c + c/) 
2(6 + c + rf) [-9a2 + 62 + (c + c/)2 - 3cd - b{c + d)] ' 
d — c 

^ " 2 [-9a2 + 62 + (c + c/)2 - 3cd - 6(c + rf)] ' ^ '^^^ 

where the overall factor y'^v'^ in these expressions has been omitted. From the discussion 
of the vacuum alignment, we know that 6 is a real parameter, and a can be real or purely 
imaginary. For the doublet flavon (p, the VEVs vi and V2 can be both real for the first 
solution in Eq. (4.15), yielding real c and d, or, for the second solution shown in Eq. (4.16), 
Vi and V2 are complex conjugates of each other, yielding c = d*. Therefore there are only 
four possible cases allowed in our model, which are listed in Table 3. 

In order to understand the phenomenological implications of the model, it is useful to 
know how the generalised CP symmetry is spontaneously broken and what the remnant 
CP symmetry is. Once the flavon fields acquire their VEVs, the 5*4 family symmetry and 
the generalised CP symmetry are spontaneously broken. Imagining that the flavon VEVs 
could transform in the way as the flavon fields under the action of CP transformation, i.e. 

('\ c. / M 

{fs)=\ 1 \vs — > p3ig){vs)* = p3'{g) \ i \v*s, g e , 

w = J ^ P2{g){<pr = P2{g) ;i j , (o = ^ ^ io* = n* , (4.26) 

then all the initial 24 consistent CP symmetries would be kept. However, the flavon VEVs 
are only numbers and they don't change at all under a CP transformation. Therefore only 
those CP transformations which transform the corresponding flavon VEVs into themselves 
remain symmetries of the theory after symmetry breaking. That means, the residual CP 
symmetry Pr{g) in the neutrino sector satisfies 

Ps\g){VsY = {Vs), P2{g)W = {4>), p^{g){0* = {0* = {0 ■ (4.27) 
In the following, we discuss the four cases in Table 3 one by one. 
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(i) In this case, both the triplet VEV vs and the doublet VEVs Vi and V2 are real. 
We can easily check that the generalised CP symmetries for h = 1 and h = S are 
preserved in the neutrino sector, i.e. the residual CP symmetry in the neutrino sector 
is H^p = {pr(l), Pr{S)}. From Eq. (4.25) we see that all four parameters a, (3, 7 and 
e are real. This is exactly the case (i) discussed in the general analysis of section 3 
which is solely based on symmetry arguments. As a result, the PMNS matrix and 
the resulting lepton mixing parameters are of the form shown in Eqs. (3.21-3.23). 
In the present case, the generalised CP symmetry Xj. = Pr(l) is preserved in both 
the neutrino and the charged lepton sector. According to the general results for 
weak basis invariants stated in Appendix A, the CP phases would be trivial, i.e. 
although this case is viable, there is no CP violation, as has been shown already in 
Eq. (3.22). 

(ii) In this case, vs is purely imaginary and the remaining VEVs vi and V2 are real. 
The generalised CP symmetry is broken to H^p = {pj.{TST'^), p^iT'^ST)} in the 
neutrino sector. Concerning the light neutrino mass matrix, the parameter a is 
imaginary, while /3, 7 and e are all real. This is exactly case (ii) of the general 
analysis of section 3, where the light neutrino masses are degenerate and hence this 
case is not viable. 

(iii) This case corresponds to the VEV vs being real together with the solution vi = V2 
for the doublet (p. Now only two of the 24 generalised CP symmetries are preserved 
in the neutrino sector and Hj^p = {pj.{U), pr{SU)}. The light neutrino mass ma- 
trix is of the form given in Eq. (4.24), where a, /3 and 7 are real while e is purely 
imaginary. This corresponds to case (iii) studied in section 3, and thus the pre- 
dictions for the PMNS matrix in Eq. (3.27) with lepton mixing parameters as in 
Eqs. (3.28,3.29). The neutrino mixing is of trimaximal form, and it is remarkable 
that we obtain maximal Dirac CP violation 6c p = ±7r/2 in this case. As for the 
generalised CP symmetry breaking, although the CP symmetries H^p = {pr(l)} 
and H^p = {pj.{U), pj.{SU)} are preserved in the charged lepton and the neutrino 
sector, respectively, the CP symmetry is completely broken in the full theory, and 
the mismatch between H^,p and H^p is precisely the origin of the maximal Dirac 
CP violation. Therefore this case is viable and predicts maximal CP violation. 

(iv) This case corresponds to vs being purely imaginary combined with the solution 
Vi = The remnant CP symmetry in the neutrino sector becomes H^p = 
{priTST'^U), pr{T^STU)}. The parameters a and e of the corresponding light 
neutrino mass matrix are imaginary, while /3 and 7 are real. This is exactly the case 
(iv) investigated in the general analysis of section 3. As has been pointed out, in 
order to achieve agreement with the present experimental data, one has to permute 
the rows of Upmns as done in Eq. (3.34). This permutation corresponds to exchang- 
ing the three charged lepton masses. However, in the present model, the charged 
lepton masses are predicted to be of different orders in the expansion parameter 
fr/A (f^/A), so that this permutation is forbidden. Although we can permute its 
columns, as the neutrino mass spectrum can be either normal or inverted, the re- 
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suiting PMNS matrix always leads to sin^ ^13 = 1/3, which is much larger than the 
measured value. Hence in the framework of the effective model, this case is not 
viable. 

Finally, we note that the residual family symmetry in the neutrino sector is Gy = 
Z2 = {1,5*} for all the four cases discussed. This is the reason why the second column 
of the tri-bimaximal mixing matrix is kept. In summary, in the effective model there 
are only two viable cases, namely (i) with no CP violation, and (iii) with maximal CP 
violation 6c p = ±7r/2. 

5 A renormalisable 5*4 >^ Hcp model 

In the effective model, the LO predictions for the flavour mixing parameters will generally 
be corrected by higher dimensional operators allowed by the imposed symmetries. As a 
result, the purely effective formulation will leave room for different physical predictions. 
Moreover, it is generally believed that the non-renormalisable operators in the effective 
theory should comes from a more fundamental underlying renormalisable theory by inte- 
grating out the heavy degrees of freedom. Therefore, the improved renormalisable model 
of this section addresses these questions while giving rise to tri-bimaximal neutrino mixing 
at leading order which is broken to trimaximal mixing by next to leading order (NLO) 
contributions. This will naturally explain why the reactor angle and the deviations from 
maximal atmospheric mixing are relatively small. 

The model of this section is inspired by the renormalisable trimaximal S4 model of [22] 
which was originally proposed without Hcp- Here we shall construct an analogous renor- 
malisable 5*4 model but with generalised CP symmetry based on 

S4 X Hcp , (5.1) 

supplemented by the same extra symmetries as in the previous effective model, i.e. 

Z4 X Z3 X U{1)r , (5.2) 

in order to control the allowed operators in the model. The matter fields, fiavon fields and 
their transformation properties under the imposed symmetries are presented in Table 4, 
while the driving fields, messenger fields, and their transformation properties are given in 
Table 5. Tables 4 and 5 may be compared to Table 2 of the effective model. 

As before we impose a generalised CP symmetry Hcp consistently with S4. In addition 
to forcing the couplings in the superpotential to be real, the generalised CP symmetry 
will be spontaneously broken. The remaining symmetry in the neutrino sector, 

G'^p = Z2 X H^p, (5.3) 

will result in predictions for 6c p and other relations, as discussed in the previous section. 
However in order to justify this breaking, we must perform a detailed analysis of the 
fiavon potential of the renormalisable model. 
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Table 4: The transformation properties of the matter fields, Higgs and flavon fields in the renormalisable 
model. 



5.1 Vacuum alignment 

The renormalisable driving superpotential that is linear in the driving fields and invariant 
under the fiavour symmetry is 

Wd = ^1 (<^r(<^T<^T)3/)i +^2 (<^r(^^T)3')l +^3C°(<^T<^t)i +^4C°(w)l 

+he (^5^5)1 + M' + fiee + m^a' + /sA^a^ , (5.4) 

where the term A^HuHd has been neglected, as it will play no role in the fiavon vacuum 
alignment because the breaking of 5*4 flavour symmetry is typically assumed to occur 
around the GUT scale. Furthermore, the couplings Qi, fi and M are real parameters due 
to the generalised CP symmetry. 

Similar to the effective model of section 4, the vacuum alignment for the flavons 
associated with the charged lepton sector is determined by the F-term conditions of Lpj- 
and C°: 



dwd 

dwd _ r, r 2 
dwd 
dwd 



This set of equations is satisfied by two solutions, where the ambiguity caused by 
symmetry transformations is ignored. The first solution is given by 




(^) = J ^ Vr^, with 4 = ' (^•^) 



(^) = J ^ Vr,, with Vt = ^Vr, . (5.6) 
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In this work, we choose the second solution since the first solution can be removed by in- 
troducing another driving field that transforms identically with similar to the effective 
model. Notice that the phase of can be absorbed into the lepton fields. Therefore, we 
take to be real, implying that the VEV vt is real as well. 

We continue the analysis of the vacuum alignment by considering the fiavon fields 
associated with the neutrino sector, i.e. ips, 4> and ^. The F-term conditions determining 
their alignments are 

dwd 
dwd 
dwd 
dwd 

dwd 

% 
dwd 

dwd 

This set of equations leads to vacuum alignments given by 



(^ll -VS2¥^S3) 




^1^82 


+ 02<^53' 


) + f3^V8, 


= , 




+ /2(<? 




+ (1^2^82. 


) + f3^V8s 


= , 




+ /2(<? 


^1^8^ 


+ 02V55i^ 


) + f3^^82 


= , 




) = o. 










h (^iV'Si - ^^32] 


) = o. 












) = o. 












+ 2/6^ 


61 02 - 


f he = 


. 





(5.7) 




where 



vl = -^{flh + Vlh)u\ v^ = -^u. (5.9) 

0/2 /5 ^/2 

Notice that the vacuum alignments of the fiavons shown in Eq. (5.8) are invariant under 
the action of both the 5* and U elements of 5*4, preserving the tri-bimaximal Klein four 
subgroup. Furthermore, the phase of u can be absorbed by field redefinition. Hence, u can 
be taken to be real without loss of generality. This renders real as well, but the VEV 
Vs can be real or purely imaginary depending on the coefficient — (/l/e + 2/I/7) /(/l/s) 
being positive or negative, respectively. 

We conclude the discussion of the vacuum alignment of the renormalisable model by 
considering the last two operators of Eq. (5.4) which are responsible for the alignment of 
the fiavon A. They provide the F-term minimisation condition 

|g = M^ + /sA^ = 0, (5.10) 
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Table 5: The driving fields, messenger fields and their transformation rules under the 84^ x Z4 x Z3 and 
U{1)r symmetries. 



which is satisfied when 

vl = -MV/8 , (5.11) 

where va = (A). Notice that if fs > 0, then va is purely imaginary. However if /§ < 0, 
then f A is completely real. Hence, f a is constrained to be either real or purely imaginary 
in this model. Thus, having completed the discussion of the vacuum alignment of the 
renormalisable model, we now proceed to investigate the leptonic masses and mixings 
predicted by it. 



5.2 The structure of the model 

The charged lepton sector is formulated at the renormalisable level with the introduction 
of three pairs of messengers Qi and Q'[ {i = 1, 2, 3). Note that these messengers are chiral 
superfields with non-vanishing hypercharge -|-2(— 2) for {^i)- With the particles and 
their transformation properties listed in Tables 4 and 5, we obtain the renormalisable 
superpotential for the charged leptons, 

+Mn, {^2^^2)1 + ^sA (^^2^^2)1' + (^^3^^3)1 , (5.12) 

where general CP invariance again implies that all the order one coupling constants Zi 
and the messenger masses M^^, M^^ and Mq^ are real. Furthermore, since the terms 
A {QiQD^, and A (^2^2)-^/ lead to corrections to the Qi and fl2 niasses, respectively, and 
the mass scales of the messenger fields are much larger than the VEVs of the fiavons, the 
contributions of these two operators can be safely neglected. 

Integrating out the messenger pairs Qi and (the corresponding Feynman diagrams 
are shown in Fig. 1), yields the following effective superpotential for the charged lepton 
masses: 

eff Z1Z2 , \ TT c , ^l^S^A ffT \ \ c 

Then, by applying the VEVs in Eq. (5.6), a diagonal charged lepton mass matrix is 
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I 
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_1_ 



I 
I 
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_1_ 



nr. 



Hi 
I 
I 
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_1_ 



I 
I 
I 

_L 



Hd 1] ri 
I I I I 
I I I I 
I I I I 
L I ^ I ^ I ^ I 

Figure 1: The diagrams which generate the effective operators for the charged lepton masses, where 
crosses indicate the mass insertions for fermions. 



obtained with 

Vt VTVr, VtV^ 

f^T = —ZiZo— — Vd-, rrir, = ZiZ-xZa— — -^^ — fw, me = —Z\Z-iZ^Ze,— — , ^ ' ^ — va. (5.13) 

The mass hierarchies among the charged leptons are reproduced naturally without invok- 
ing another mechanisms. Here the VEVs of the flavons ^px and r] are responsible for the 
spontaneous breaking of both flavour and generalised CP symmetries. As the effective 
model in section 4, the 5*4 flavour symmetry is broken completely and the remnant CP 
symmetry is H^p = Pr(l) in the charged lepton sector. 

Having completed the analysis of the charged lepton sector, we now turn to discuss the 
neutrino sector. We begin this by writing the renormalisable superpotential responsible 
for the light neutrino masses, which consists of the LO and relevant messenger terms: 



w, = w^'' + w^, (5.14) 



where 



= xi((iV'=S)3,(^5)i + X2((iV'=S)2</.)^+X3(iV^S^)i, A + Ms(SS^)i . (5.15) 

where the messenger field S (S'^) is a chiral superfield carrying zero hypercharge, and all 
the parameters Xj and yi are real due to generalised CP invariance. It is clear that the 
Dirac neutrino mass matrix takes the simple form 

1 

mn = yvu \ 1 \ . (5.16) 
1 
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I I I 

I I I 

I I I 

Figure 2: The diagrams for the RH neutrino masses, where the cross indicates a fermionic mass insertion. 



Further notice that w^^ gives rise to the RH Majorana neutrino mass matrix m^^, and 
after inserting the flavon VEVs of ips-, <p and ^ from Eq. (5.8), this is revealed to be 



^M=yi^. I -1 2 -1 \+y2u\ 1 \+yzvA 1 1 1. (5.17) 
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The resuhing effective hght neutrino mass matrix m^*^ = —m£,{m^f) is exactly 

diagonalised by the tri-bimaximal mixing matrix Utbi 

U^bK^'Utb = diag (mf«, m^^, m^«) , (5.18) 

where the light neutrino masses rn^2 3 



1^1 = -T, \ , ^2 = ' ^3 



^yivs + y2U - ysv^ y2U + 2y3v^ 3yiVs - y2U + y^v^j, 

(5.19) 

The reason why the tri-bimaximal mixing is produced is because the VEVs oi^ps-, (p ^"^^ ^ 
preserve the Klein four subgroup generated by the tri-bimaximal S and U generators, as 
has been pointed out in section 5.1. The LO and NLO contributions to the RH neutrino 
masses are shown diagrammatically in Fig. 2. Integrating out the messenger fields S and 
S*^ yields the NLO effective operator 

<"^ = -^A((iV'^iV=), </>),, . 

The effective operator A {{N'^N'^)^ips)ii is not reproduced, since the contraction N'^N'^ 
vanishes from the antisymmetry of the associated 5*4 Clebsch-Gordan coefficients, as is 
shown in Appendix B. Straightforwardly, we see that w^^^ gives rise to a NLO contri- 
bution to niM of the form 

NLO V^V^ 

Note that the NLO contribution mj^^'^ is induced by the flavon A which further breaks the 
remnant Klein four symmetry to Then, the RH neutrino mass matrix rriM including 
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the NLO contribution can be expressed as 

TUm — 




(5.20) 



where the parameters a, b, c and d are defined as 



a = yivs, b = y2U, c = y^v^, d = X2X^-^ . (5.21) 

My, 

The hght neutrino mass matrix mi, is given by the seesaw formula 




It is the most general neutrino mass matrix consistent with the residual fiavour sym- 
metry, as is shown in Eq. (3.12). The parameters a, /3, 7 and e are given by, 

1 2(6 - 5) 



-9S2 + (6 - 5)2 + 3J2 ' ^ 3(6 + 25) '''3 



9a2 - (6 - 5)2 _ 3J2 
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1 b-c d 



3(6 + 25) 3 9^2 _ (5 _ 5)2 _ 3J2 _9a2 + (6 -5)2 + 3^2 



(5.23) 



where the overall factor y'^v'^ has been omitted. We note that the first three terms in 
the light neutrino matrix of Eq. (5.22) preserve tri-bimaximal mixing while the last term, 
proportional to e, violates it. In the present model, the e term is induced by the NLO 
contributions and hence suppressed by va/My with respect to a, /3 and 7. This provides 
a natural explanation as to why the reactor angle as well as the deviations from maxi- 
mal atmospheric mixing are relatively small. Although, their definite values cannot be 
predicted. 

From the vacuum alignment in section 5.1, we know that both b and 5 are real pa- 
rameters, and a and d can be real or purely imaginary. Therefore, there are four possible 
cases allowed by the domains of the parameters in the renormalisable model, these are 
summarised in Table 6. It turns out that each of these four cases will preserve different 
generalised CP transformations as well as lead to different phenomenological predictions. 
In what follows, we shall discuss the four cases in Table 6, one by one. 

(i) In this case, both the triplet VEV vs and the singlet VEV are real. The gener- 
alised CP symmetries for g = 1 and g = S are preserved in the neutrino sector, i.e. 
the residual CP symmetry in the neutrino sector is H^p = {pr{l) , pr{S)} , where r 
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a 


/3 


7 


e 


(i) 


aeR, 6eR, £eR, JeM 


real 


real 


real 


real 


(ii) 


aeC, 5eM, ceK, deR with a = -a* 


imaginary 


real 


real 


real 


(iii) 


aeR, cgR, dec with d = -d* 


real 


real 


real 


imaginary 


(iv) 


aeC, feeR, ceR, JeC with a = -d* and d = -d* 


imaginary 


real 


real 


imaginary 



Table 6: The four allowed ranges of a, /?, 7 and e as dictated by the domains of a, b, c and d in the 
renormalisable model. 



denotes the irreducible representations of 5*4. For the corresponding hght neutrino 
mass matrix shown in Eq. (5.22), the parameters a, /3, 7 and e are all real, yielding 
a matrix that is precisely the same as the neutrino mass matrix of case (i) from the 
general analysis in section 3. Therefore, the neutrino mixing matrix is of trimax- 
imal form, and the predictions for lepton mixing angles and CP phases are given 
by Eqs. (3.22,3.23). There is no CP violation, as the CP symmetry = Pr(l) is 
conserved in both the neutrino and charged lepton sectors. However this case gives 
a viable description of the lepton mixing angles. 

(ii) In this case, vs is purely imaginary and va is real. The generalised CP symmetry 
is broken to H^p = {pj-iTST"^), pr{T^ST)} in the neutrino sector. The resulting 
parameter a is imaginary and f3, 7 and e are real. As a result, this is case (ii) 
discussed in section 3, and the light neutrino masses are degenerate, i.e. \mi\ = Im^l. 
Hence, this case is not viable. 

(iii) This case corresponds to the VEV vs being real and being purely imaginary. 
Only two of the 24 generalised CP symmetries are preserved in the neutrino sec- 
tor, i.e. H^p = {pr{U), pr{SU)}. Regarding the light neutrino mass matrix of 
Eq. (5.22), the tri-bimaximal violating parameter e is imaginary and a, /3 and 7 are 
real. Therefore, this case is identical to case (iii) of the general analysis inspired 
by symmetry arguments, and the predictions for its mixing angles and CP phases 
are as given in Eqs. (3.28,3.29). Notice that this case produces maximal Dirac CP 
violation \6cp\ = 7r/2, which is the result of the mismatch between the remnant CP 
symmetries H^p = {pr(l)} and H^p = {pr{U), pr{SU)}. This case is therefore 
viable. 

(iv) Both the VEVs vs and f a are purely imaginary in this case, the remnant CP sym- 
metry in the neutrino sector turns out to be H^p = {pj.{T ST'^U) , p^iT"^ STU)} . The 
parameters a and e are imaginary while /3 and 7 are real. Hence, the corresponding 
light neutrino mass matrix has the same form as that of case (iv) analysed in section 
3, and the resulting PMNS matrix is given in Eq. (3.33) up to permutations of rows 
and columns. However, in the present model, the mass hierarchies among the three 
charged lepton are reproduced, as is shown in Eq. (5.13). Consequently, permuting 
the rows of the PMNS matrix is forbidden. Notice that even if we desired to ex- 
change its columns, as the neutrino mass order is less constrained so far, we would 
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always get sin^ ^13 = 1/3, which is much larger than the experimental observations. 
Therefore, we conclude that this case is not viable for the renormalisable model. 

In summary, in the renormalisable model there are only two viable cases, namely (i) 
with no CP violation, and (iii) with maximal CP violation. This is the same as in the 
effective model. 

6 Conclusions 

The measurement of the reactor mixing angle, which is observed to be rather large, offers 
some encouragement that the measurement of leptonic CP violation, in particular the 
Dirac oscillation phase 5cPi niay be possible in the not too distant future. This has led 
to renewed interest in theories that are able to predict the value of 5c p- 

In this work, we have focused on a promising framework where (generalised) CP sym- 
metry in the lepton sector is combined with a discrete family symmetry (and perhaps 
other symmetries). The CP and family symmetries are subsequently spontaneously bro- 
ken, with different subgroups preserved in the neutrino and charged lepton sectors. The 
models presented in this paper provide a generalisation for the direct models reviewed 
in [8] to include CP symmetry. 

Within this framework, we have addressed the question of spontaneous breaking of a 
generalised CP symmetry in particular for models based on an 5*4 family symmetry. We 
have constructed two models of leptons based on family symmetry combined with a 
generalised CP symmetry Hcp-, one at the effective level and another one at the renormal- 
isable level. In both models, we have shown how the flavon potential can spontaneously 
break the symmetry 5*4 x Hcp down to Z2 x H^p in the neutrino sector. This symmetry 
breaking was simply assumed to happen in [20] without any dynamical justification. 

In our models, the choice of preserved CP symmetry H^p in the neutrino sector is 
controlled by free (real) parameters in the flavon potential, enabling us to dial the type of 
CP violation. Of the two realistic models of this kind that we have proposed, one at the 
effective level and another one at the renormalisable level, we find that both models predict 
trimaximal lepton mixing with CP being either fully preserved or maximally broken and 
the intermediate possibility forbidden by the structure of the models. 
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Appendices 



A Implications of a generalised CP symmetry 

It is well-known that at low energies there are three CP violating phases in the lepton 
sector: one of the Dirac type and two of the Majorana type. The strength of leptonic 
CP violation of the Dirac type, which can be observable through neutrino oscillations, is 
determined by the following CP-odd weak basis (WB) invariant [24]: 

J'l = Tt [hi, hi]^ , with = niym\ and hi = raim\. (A.l) 

It can be fully written in terms of physical observables as 

J^ = -Qi {ml - ml) {ml - ml) {ml - ml)^m\^ ^m\ Am^^ ^cp, (A.2) 

where Am^^- = — are the light neutrino mass squared differences and the quantity 
J CP is the Jarlskog invariant [25], 

JCP = Im[{UpMNs)n{UpMNs)22iUpMNs)l2iUpMNs)2l\ 

= - sin(26'23) sin(26'i2) sin(26'i3) COS 6113 sin ^cp- (A.3) 

o 

where the PMNS matrix Upmns has been parameterised as in Eq. (3.17), 6ij are the 
mixing angles chosen to lie in the first quadrant, and Sep is the Dirac CP-violating phase. 
It is obvious that J7i vanishes for 6cp = 0, and vice versa the Dirac phase 6cp would be 
zero if J7i vanishes. From Eq. (2.3), we see that the generalised CP invariance implies 

hi = XLhfxl K = XlhlXl. (A.4) 

Therefore, we have 

Ji = Tr [hi, hif = Tr [x^ [K, hjfx^^ = Tr [(/i^)^, /if]' = -Tr [hi, = 0. (A.5) 

This indicates that there is no Dirac type CP violation if the generalised CP symmetry is 
preserved. In the case of Majorana neutrinos, there is also the possibility of the Majorana 
type CP violation. It has been established that the vanishing of the WB invariant Ji 
together with the following two WB invariants, 

J2 = Im Tr {hihlmlhlmy) , 

J3 = Tr [mi^himl, h*if , (A. 6) 

provides necessary and sufficient conditions for low energy CP invariance [21,26-28]. No- 
tice that and jTs measure the Majorana type CP violation. Then, from the generalised 
CP invariance requirement of Eq. (2.3), we have 

J2 = Im Tr (^XLhlh^m^himlXlj = Im Tr {hihlmlh*m^)* = 0, (A.7) 

J3 = Tr {XI [mlh]m„ hif Xl) = Tr [{m^ml)'^ , {hl)^f = -Tr [m^ml, hlf = 0. 

This means that CP will be conserved if the theory is invariant under a generalised CP 
transformation. As a result, the generalised CP symmetry must be spontaneously or 
explicitly broken in order to generate non-trivial CP phases. 
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1, 1' 



s 



T 



U 
±1 



1 
1 



uj 
u;2 



1 

1 



3, 3' 



-1 





Table 7: The iS'4 representation matrices for the S, T and U elements in different irreducible representa- 
tions, where u) = e^'^'/'^. 



B Group theory of S'4 

5*4 is the permutation group of four distinct objects, and it is isomorphic to the symmetry 
group of a regular octahedron. It has 4! = 24 elements and can be expressed in terms of 
two generators. However, in order to clearly see the connection to the groups ^4 and 
it is useful to express S'4 in terms of three generators, 5, T and U [29], which obey the 
multiplication rules 

= T^ = U^ = {STf = {SUf = {TUf = {STUY = 1 . (B.l) 

Notice that the generators S and T alone generate the group ^4, while the generators 
T and U alone generate the group S3. Taking all possible combinations of S, T, and 
U (subject to the rules of Eq. (B.l)), yields the 24 elements of S'4 which belong to 5 
disjoint conjugacy classes. To emphasise the geometric aspect of S'4, we adopt Schoenflies 
notation, in which kCn designates a conjugacy class of k elements that are all rotations 
by to express the conjugacy classes of S'4 as 

ICi = {1}, 

3C2 = {S,TST^,T^ST} , 

QC'^ = {U, TU, SU, T^U, STSU, ST^SU} , (B.2) 

8C3 = {T,ST,TS,STS,T^,ST^,T^S,ST^S} , 

6C4 = {STU, TSU, T^SU, ST'^U, TST'^U, T'^STU} . 

These conjugacy classes can be used to deduce the various irreducible representations of 
S'4 because (by theorem) the number of the irreducible representations must equal to the 
number of conjugacy classes. Thus, S'4 has five irreducible presentations. Additionally, 
the sum of the squares of the dimensions of the irreducible representations must equal the 
order of the group, i.e. 24. This implies the five irreducible representations of S'4 are two 
1-dimensional (1, 1'), one 2-dimensional (2), and two 3-dimensional (3 and 3') irreducible 
representations. Our choice of the explicit basis for the representation matrices of S, T 
and U is listed in Table 7 [22]. Furthermore, the Kronecker products of these 5 different 
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irreducible representations are 

1®R = R, 1'®1' = 1, l'®2 = 
2 ® 2 = (1 © 2)s © i; , 2 ® 3 = 2 ® 
3©3 = 3'©3' = (1©2© 3')s © 3a , 



2 , 1' © 3 = 3' , 1' © 3' 
3' = 3 © 3' , 
3 © 3' = 1' © 2 © 3 © 3' , 



3, 



(B.3) 



where R stands for any 5*4 representation, and the index s (a) denotes symmetric (anti- 
symmetric) combinations. These Kronecker products, along with the explicit forms of the 
generators in Table 7 can be used to calculate the corresponding Clebsch-Gordan (CG) 
coefficients. These CG coefficients can be also be found in other works, e.g. Ref. [22], but 
we list them here for completeness. In the following reporting of the CG coefficients of S4, 
we use Oi to denote the elements of the first representation and f3j to indicate those of the 
second representation of the product. Furthermore, "n" counts the number of "primes" 
in the Kronecker product (e.g. in 1 © 1' = 1', n = 2). 
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3®3 ^3'] /2ai/3i -02/33-03/32 
3(0 (g, 3(') ^ 3(0 <( n = odd 3 ® 3' ^ 3 > 2^3/33 - ai/32 - a2/3i 

3'® 3' ^ 3' ) \2a2/32 - a3/3i - ai/33^ 

I - ^ - / 02/33 - "3/32 

3(') ^ 3(0 ^ 3(0 <j n = even 3' ® 3' ^ 3 J- ai/32 - a2/3i 

\a3/3i - ai/33^ 

Now that we have discussed the simpler aspects of the group theory of 5*4, we turn to 
the more complex topic of the automorphisms of the group 6*4, i.e. the generalised CP 
transformations X,.. 

Although the most general solution for the Xj. is given by the representation pr{h) 
of the 5*4 group elements h G 5*4, as is claimed in Ref. [17], it is still instructive to see 
how this explicitly occurs. Let us first consider the consistency equation Eq. (2.7) for the 
generators S, T and U in the faithful representations 3 and 3': 

^3{3')P3(3')('^)^3(3') = P3(3')('S"), 

-^3(3')P3(3')(^)-^3{3') = P3i3'){T'), (B.4) 
^3(3')P3(3')(^)^3(3') = P3{3'){U'). 

It is clear from the above equations that the orders of primed and unprimed generators 
must be identical, i.e. 5" and U' must be order 2 elements and T' an order 3 element. 
As a result, the consistency equations constrain the mappings (automorphisms) such that 
the generators S, T, U can only be mapped into specific (unions of) conjugacy classes. 
Namely, 

S' e 3C2 U , T' G 8C3 , U' e 3C2 u QC'2 . (B.5) 

Additional constraints can be derived by considering the singlet representation 1'. Since 
S = T = 1 and U = —1, S' cannot be in the conjugacy class GCg. Likewise the element 
U' cannot be in the conjugacy class 3C2. This leaves us the automorphisms 

S' G 3C2 , T' G 8C3 , U' G . (B.6) 

By investigating all 3 x 8 x 6 = 144 possible values of 5", T' and U', we find only 24 
solutions that satisfy the corresponding consistency equations Eq. (B.4), which can be 
compactly written as 

^3(3') = p3{3'){h), he S4 , (B.7) 

with 

^3(3')P3(3')('^)^3(3') = P3i3'){hSh~^), 
^3(3')P3(3')(^)^3(3') = P3{3'){hT'^h~^), 

^3(3')P;(3')(^)^3"(3') = P3i3'){hUh~'). (B.8) 

Continuing to the 2-dimensional irreducible representation, we have 

PliS) = P2{S), pliT) = p2{T'). pUU) = P2{U). (B.9) 
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Thus, the generahsed CP transformations consistent with the faithful three dimensional 
irreducible representations are 



X2 = Piih), h e S4 . 
Finally, for the singlet representations 1 and 1', we take 

^i(i') = Pi(i')(^)' h e S4 . 



(B.IO) 



(B.ll) 



Therefore the generalised CP transformation consistent with an S4 family symmetry is of 
the same form as the flavour group transformation, i.e. 



Xr. 



Pr{h) 



he Si 



(B.12) 



This provides a proof to the statement in Ref. [17] that the generalised CP transformation 
group consistent with a 5*4 flavour group is 5*4. 

We end this appendix by noting that our basis for S, T and U is related to the basis 
choice in Ref. [20] (denoted here as S*, T and U) by the unitary transformation 



V 

such that (for the 3' representation) 
S = V^SV = 



\ 



2 J_ 

3 V3 

1 1 i_ 

\ L J_ ^ 

\ Ve Vs V2 / 



(B.13) 





T = V^TV = - 
2 



U = V^UV 

and the generalised CP transformations in the two bases are related by 

-^^3(3') 

It is interesting to note that the matrix V satisfies 





(B.14) 



(B.15) 




(B.16) 



so that ^^3(3/) = V''Xs(s')VU is explicitly an element of S4 in the basis with the tilde, 
given that Xs{s') is an element of S'4 in our basis. 



^Said again, the generalised CP transformation group is the identity up to inner automorphisms. 
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